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Abstract. - We demonstrate the possibility of existence of meta-stable JV-body Efimov states 
in trapped Bose systems with large scattering length. We calculate spectra of trapped systems 
of iV=3,4,5,6, and 7 bosons using a stochastic variational method with a restricted correlated 
Gaussian basis. For each system the calculations reveal a series of Efimov states where the energy 
and the r.m.s. radius exhibit the characteristic exponential dependence upon the state number. 
We also estimate the contribution of these states to the recombination rate of Bose-Einstein 
condensates. 



£>y Introduction. — The Efimov effect [1] appears in 
, <"*] quantum three-body systems when attractive interactions 
Qjjetween at least two pairs of particles are such that the 
scattering length is much larger than the range of the in- 
y—\ teraction; in other words two of the three two-body sub- 
!>■ systems are close to the threshold of binding. Under these 
conditions a characteristic series of weakly bound and spa- 
^^tially extended states, called Efimov states [2], appears 
C*~)in the system. These states appear due to specific long 
range two-body correlations between particles caused by 
j^^the large scattering length. 

qq The effect is easiest to see in the hyper-spherical adia- 
^^batic approximation where the slow adiabatic variable is 
^the hyper-radius p (root square radius of the system). It 
has been shown [2] that close to the two-body threshold 
the effective adiabatic potential W(p) is attractive 1 and 
^asymptotically proportional to the inverse square of the 
" hyper-radius, 
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where m is the mass scale, and £ is a constant depending 
on masses of the particles [3,6]. 

A sufficiently large positive £ 2 leads to a geometric series 
of bound states with exceedingly small energies, E n oc 
e~^' n , and exceedingly large root mean square radii, R n oc 
ea^'™, where n is the state number and £ = 

There is a general theoretical consensus [4-6] that the 
first excited state of the helium trimer 4 He3 is an Efimov 
state, although the experimental observation so far proved 
elusive [7]. 

The three-body Efimov effect has recently got some sup- 



port [8] from an experiment with trapped Bose gases where 
the recombination rate was measured as function of the 
scattering length. The latter was varied using the Fes- 
chbach resonance technique by applying an external mag- 
netic field. A sharp peak in the recombination rate was 
detected and interpreted as a three-body Efimov resonance 
in qualitative agreement with theoretical predictions [9]. 

It has been shown that in an TV-body system with N >3 
the Efimov effect docs not exist at the A— 1 threshold [10]. 
At the two-body threshold the A-body Efimov states with 
N > 3 can not exist either as the clusters with 3 and 
higher number of particles are generally deeply bound at 
this point. The Efimov states would then be unstable 
due to lower lying thresholds and would decay into deeply 
bound cluster states. 

However it has been suggested in [11] that a sequence 
of meta-stable iV-body states with the characteristic ex- 
ponential energy dependence can yet show up at the two- 
body threshold. Using the TV-body hyper-spheric method 
it has been shown that an iV-body system at the two- 
body threshold has a hyper-spheric adiabatic potential 
with inverse-square dependence. This peculiar adiabatic 
potential appears due to the same mechanism as for three 
particles and thus gives rise to iV-body Efimov states with 
a structure similar to that of three-body Efimov states: an 
(otherwise) uncorrelated system with very specific two- 
body correlations caused by the large scattering lengths. 

This specific hyper-spheric adiabatic potential is not the 
lowest one as different bound clusters with lower thresh- 
olds create lower lying adiabatic potentials. However, 
although not truly bound, these ./V-body Efimov states 
might still exist as meta-stable states slowly decaying into 
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clusters, much like the Bose-Einstein condensate states. 
The structure of the Efimov states is determined by the 
long-range two-body correlations and is thus quite dissim- 
ilar to the structure of clusterized states with short-range 
many-body correlations. Therefore the overlap between 
Efimov states and clusterized states should be small and 
consequently the life-time should be large. 

The conclusions about meta-stable TV-body Efimov 
states were obtained in [11] in an extreme hyper-spheric 
adiabatic approximation where the couplings to all other 
channels were neglected. In this letter we report on a more 
realistic calculation of A-body Efimov states with a dif- 
ferent method, namely the restricted correlated Gaussian 
stochastic variational method where no adiabatic approx- 
imation is assumed. 

The system and the method. — We consider a sys- 
tem of A identical bosons with mass m and coordinates 
in a spherical harmonic trap with frequency lo, where the 
scattering length is assumed to be tuned to a large value, 
facilitating the Efimov effect. 

The Hamiltonian of the system is given by 
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where the system parameters are m = 86.909 u, and the 
trap length b t = \Jfij (tow)=46189 au. 

The two-body potential is an attractive Gaussian with 
the range fo=11.65au and depth 1.24825au corresponding 
to the scattering length of a=— 1.4 • 10 6 au 3> b t . 

The wave-function of the system is represented as a lin- 
ear combination of K basis-functions taken in the form of 
symmetrized correlated Gaussians, 
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where the total angular momentum is zero, S is the sym- 
metrization operator, and Ck and ot^ are variational pa- 
rameters. The linear parameters Ck are determined by 
diagonalization of the Hamiltonian while the non-linear 
parameters af^ are optimized stochastically [14, 15] by 
random sampling from a region that covers the spatial 
distances from b to bt- 

During calculations of a given system the number of 
Gaussians in the basis is increased and the stochastic op- 
timization is carried out until all energy levels of interest 
are converged. 

With uncorrelated Gaussians the method is equiva- 
lent [15] to a mean-field approximation which is unable 
to describe the Efimov effect and which has the validity 
condition na 3 <C 1. Introducing successive correlations be- 
yond the mean-field improves this validity condition and 
allows calculations of the Efimov states. 



The structure of the A-body Efimov states is analogous 
to that of the three-body Efimov states: the spatial ex- 
tension of these states is much larger than the range of 
the potential 6, hence the density of the system is small, 
nb 3 <C 1; and there are no cluster substructures. For low- 
density systems, nb 3 <C 1, only two-body correlations are 
of importance. 

Allowing only two-body correlations, the variational 
wave-function can be simplified as 
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= S^Ck 



k=l 

x eX p(-la(V-5/3 (fe) (ri-r 2 ) 2 ) , (4) 

where p 2 = J2i<j ( r i~ r j) 2 1S the hyper-radius and aS k > and 
p( k ^ are the nonlinear parameters. The symmetrization of 
this function can be done analytically [15] which greatly 
simplifies the numerical calculations. 

This form of the variational space proved very success- 
ful in describing the Bose-Einstein condensate states [12] 
which have similar structure. The clusterized states, un- 
interesting in the present context, are explicitly excluded 
while the crucial two-body correlations responsible for the 
existence of the Efimov states are retained. Otherwise it 
would be impossible to calculate the Efimov states due 
to the enormous number of cluster thresholds and cluster- 
ized states in an A-body system. Omitting the clusterized 
states, into which the Efimov states would decay, effec- 
tively disregards the widths of the Efimov states. However 
the clusterized states differ significantly in spatial struc- 
ture and therefore should overlap very little with the com- 
puted Efimov states. Indeed an estimate [13] of the life- 
times of the molecular three-body Efimov states showed 
that their widths are extremely small and that they can 
be considered bound for all practical purposes. 

Results. — The energies of Efimov states in a trap 
should be on one hand much larger than the typical energy 

fc2 

scale of a cluster state, -^W, and on the other hand smaller 
than the oscillator energy Too — 

Wc have calculated the spectrum of trapped boson sys- 
tems with the Hamiltonian eq. (2) with A=3,4,5,6, and 7. 
The calculated energies are shown on Fig. 1. Indeed in 
the indicated energy region for each of the A-body sys- 
tems there is a series of states with exponential depen- 
dence upon the state number, E n cx e~^ Nn . 

The exponential fits give the numbers C3=6-33±0.032, 
C 4 =3.40±0.14, Cs=l-79±0.034, C 6 =1.31±0.020, and 
^7=1.01±0.007. The value for A=3 agrees within 3 sigma 
with the known analytical result of 6.244 (see [3] and refer- 
ences therein) . Fig. 2 shows our numerical results together 
with the asymptotic expression for large A, 
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Fig. 1: Energy E n as function of the state number n for nega- 
tive energy states of a system of N bosons in a harmonic trap 
with length bt interacting via an attractive potential of range 
b with scattering length much larger than b t . The horizontal 
lines demarcate the region where the Efimov states can exist in 
trapped gases. The lines show the exponential fits of the form 
E n oc e -1 -'™ drawn through the points in the indicated region. 



which was established in [16]. Our results seem consistent 
with the asymptotic estimate. 

The spatial extension of Efimov states in trapped sys- 
tems must be much larger than the interaction range b and 
much smaller than the trap length bt- On Fig. 3 are shown 
the calculated r.m.s. radii R n as function of state number 
n. The radii of the Efimov states identified on Fig. 1 are 
reproduced well with the exponentials R„ oc e^^" where 
the parameters Cat are taken from the fits on Fig. 1. Ap- 
parently all these states fall within the correct boundaries 
and the values of the radii follow the correct exponential 
trend. 

There are several states in the 7V=6 and 7 systems with 
radii much smaller than those of the typical states in the 
series with similar energies. Clearly these states are not 
Efimov states but rather relatively compact states with a 
different structure. 

Recombination reactions into shallow TV-body 
states. — The TV-body Efimov states could in principle 
be identified by their contribution to the recombination 
rate of a cold gas as function of scattering length simi- 
lar to the three-body case [8] . As the scattering length is 
increased using the Feschbach technique, the TV-body Efi- 
mov states crossing the threshold should produce peaks 
in the recombination rate. Since the n-th Efimov states 
appears when the scattering length is about a « aoe 5 ^™ 
(where a is the scattering length corresponding to the 
lowest Efimov state) the peaks will appear as an exponen- 
tial sequence as function of the scattering length [17]. 



Fig. 2: The exponents £jv from Fig.l as function of the boson 
number TV together with the asymptotic form eq. (5). 

The scale of the rate can be estimated as follows. Let 
us first consider a three-body reaction. The rate of the 
loss of particles from a cold gas due to the three-body re- 
combination reaction into a shallow dimcr with the energy 
is given by Fermi's golden rule, 
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where the factor 3 is there since each recombination re- 
action removes three particles from the cold gas, TVq/6 is 
the number of triples in the gas of TVo particles, dvf is the 
number of final states (dimer plus the third particle) with 
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and Tfi is the transition matrix element from the initial 
three-body to the final dimer+particle state. 

In the typical experimental regime, where the scatter- 
ing length is still much smaller than the size of the trap, 
a <C bt, the transition matrix element for the non-resonant 
three-body recombination rate from a cold gas state into a 
shallow dimer state can be estimated perturbatively sub- 
stituting the asymptotic expressions for the initial and fi- 
nal wave-functions, 



Tfi = fd s rd s R[Mr)^ 
(U(R- |r) + C/(R+±r)) 



e' kr e"" 1 



(8) 



where r is the distance between two particles, R is the dis- 
tance between their center of mass and the third particle, 
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Fig. 3: Root mean square radius R n as function of the state 
number n of a system of N bosons from Fig.l. The lines show 
the fitting curves of the form R n oc e'^™ drawn through the 
same points and with the same parameters (n as on Fig.l. 



V oc b 3 is the normalization volume, k and q (k ~ q oc 
67 1 <C a" 1 ) are the initial momenta of the cold gas parti- 
cles, i/jd(f) is the s-wave function of the shallow dimcr with 
the binding energy . Using the zero-range approxima- 
tion for the transition interaction, U(r) = 47r ^ a S(r), and 

for the dimer wave-function, ipd(r) oc a" 1 ! 2 - — — , the ma- 
trix element (8) in the limit k ~ q <C a -1 is estimated as 
(cf. [18]) 

T n oc • 0) 

Finally, the non-resonant factor of the three-body re- 
combination rate becomes (cf. [3,18]) 



might be substantially different from the shallow A-l clus- 
ter plus one particle. 

Although in the regime na 3 -C 1 the A^-body recombi- 
nation has an additional small factor (na 3 ) N ~ 3 it might 
still be possible to observe the A-body Efimov states as 
a sequence of resonant peaks in the recombination rate 
as function of scattering length. The ratio of scattering 
lengths corresponding to the adjacent peaks, say number 



1 and number n, is given as 
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Conclusion. — We have calculated the spectrum of 
trapped A-boson systems with A=3,4,5,6, and 7 using 
the stochastic variation method with a restricted corre- 
lated Gaussian basis. Only two-body correlations were 
allowed in the variational space. Thus the cluster states 
were a priori mostly excluded from the variation space, 
which made the calculations technically possible. 

For each system a series of states is found with specific 
exponential dependences of the energies and r.m.s radii 
on the state number, which is a characteristic feature of 
Efimov states. For the A=3 system the exponent obtained 
agrees well with the known analytical value and the trend 
also agrees well with the existing asymptotic estimate. 

Inclusion of the cluster states would turn the Efimov 
states into meta-stable states. However the life-time of 
these states should be comparable with the Bose-Einstein 
condensate states which have a similar structure and sim- 
ilar decay modes. 

It might be possible to observe the 4-body Efimov states 
as peaks in the recombination rate of a condensate as func- 
tion of the scattering length, similar to the 3-body case. 
However, for a dilute gas the 4-body recombination rate 
is a factor no? smaller than the 3-body rate, therefore the 
accuracy constraints on the experiment would be higher. 

In conclusion we have lent theoretical support to the 
possibility of existence of long lived meta-stable A-body 
Efimov states in trapped Bose gases. 
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where n = Nq/V is the density. 

For an A-body recombination reaction into a shallow A- 
1 body Efimov state with the binding energy of the order of 
— — t the modifications to the expression for the rate (6) in- 
elude the extra 3(A — 3) spatial dimensions in the integral 
in eq.(8), which gives an extra factor (a 3 / 2 !/ -1 / 2 ) 2 ^ -3 ) , 
and also the factor N$/6 is substituted by /NL Thus 
for the A-body recombination rate we have 
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We emphasize that these simple estimates only refer to 
the non-resonant contributions and they also do not in- 
clude other types of decays where the final state structure 
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